Nature of a topological quantum phase transition in a chiral spin liquid model 
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We study the finite temperature nature of a quantum phase transition between an Abelian and a 
non-Abelian topological phase in an exactly solvable model of a chiral spin liquid pQ. By virtue of 
the exact solvability, this model can serve as a testbed for developing better measures for describing 
topological quantum phase transitions. We characterize this phase transition in terms of the global 
flux and entanglement entropy, and discuss to what extent the existence of a topologically ordered 
ground state with non-Abelian excitations is revealed at finite temperature. 



Characterizing and detecting topological order is one 
of the central questions in the field of topological phases. 
The challenge lies in that these new type of quantum 
ground states are not associated with any local broken 
symmetry. Of broader interest in the context of quantum 
phase transitions (QPT) is a question of the nature of a 
quantum critical point when a system enters a topologi- 
cally ordered phase [2]. In a conventional QPT, in which 
a local order parameter starts to gain an expectation 
value at the quantum critical point (QCP), the nature of 
the QCP is of ultimate significance. Even though it is ef- 
fectively a point of measure zero, it governs a much larger 
phase space often called the "quantum critical region" 
(See Figjl[a) where we denote the expectation value of 
an order parameter by (</>)). One can ask if, and to what 
extent, an analogy holds for topological quantum phase 
transitions. For such a question, we need a formulation 
and understanding of measures of topological order at 
finite temperature. 

Since Wen and Niu [3] coined the term "topological or- 
der" in association with the ground state degeneracy of 
(Abelian) fractional quantum Hall (FQH) states on topo- 
logically non-trivial surfaces, such ground state degener- 
acy has been widely used as an indicator of topological 
order including the non-Abelian FQH states [I]. Further, 
the implications of such degeneracy on fractionalization 
has also been discussed [5j [B]. However, the extension of 
this indicator, which is defined at T = and not directly 
accessible experimentally, to a measure at finite temper- 
ature is an open question. 

More recently, the concept of "topological entangle- 
ment entropy" has been gaining interest as an indicator of 
topological order [7J[5] or a topological QCP The cor- 
responding quantity at finite temperature has also been 
studied [10]. However, one of the issues with topological 
entanglement entropy is that it does not always distin- 
guish phases with obviously different topological orders 
such as weak pairing (vortices follow non-Abelian statis- 
tics) and strong pairing (vortices are Abelian) p + ip 
superconductors [11]. This shows that the information 
about the ground state is significantly condensed upon 
mapping to a single entropic quantity. 

We start with the observation that a chiral spin liquid 



(CSL) model [T] can provide an ideal testbed for develop- 
ing a better understanding of a topological QPT by play- 
ing the role that the transverse field Ising chain played 
in the study of conventional QPTs. It has the virtue of 
being exactly solvable, and exhibiting a non-trivial QPT 
between non-Abelian(nA) and Abelian(A) phases anal- 
ogous to the weak pairing and strong pairing limits of 
a p + ip superconductor; the same physics can also be 
accomplished for the honeycomb model with three-spin 
interaction [13]. We take a twofold approach: First, we 
employ the notion of an "expectation value" of a global 
flux operator introduced by Nussinov and Ortiz [12] as 
a finite temperature extension of the concept of ground 
state degeneracy. Second, we contrast this result to what 
can be learned from entanglement entropy. 
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FIG. 1: Phase diagrams in g — T phase space. (a)A typical 
QPT phase diagram for conventional order with the quan- 
tum fluctuations controlled by tuning parameter g. Here we 
sketched a case with the dynamical critical exponent z = 1. 
(b)Topological QPT phase diagram. 

Model- The exactly solvable CSL model on the star 
lattice PP is a variant of a spin model with topological 
order first introduced by Kitaev[T3] on the honeycomb 
lattice. In this variation, ground states spontaneously 
break time reversal symmetry and a QPT between A and 
nA phases is accessible through the exact solution. For 
brevity we employ a Majorana fermion representation of 
the model. We represent spin- 1/2 Pauli operators erf 
(a = x,y,z) of the original spin model at each lattice 
site ip] by four species of Majorana fermions Ci and df, 
erf = icidf under the constraint 

A = cid^df = 1 (1) 



so that a x a y a z = i as is expected of spin 1/2 operators. 
In terms of these Majorana fermions, the Hamiltonian is 



H[{tr ij }] = j £ 

x,y,z — link 



UijZCiCj , (2) 



a;', y', z' — link 



where t/j. 



-idfd'j is defined at each a and a' bonds 



between sites (see Fig(2]) and acts as a Z 2 gauge 

field living on the ij bond [18] . As C/y has no dynamics 
([C/jj, = 0) it can be replaced by a set of Z 2 variables 
Uij = ±1 reducing Eq. ([2| to a quadratic Hamiltonian 
7i[{ujj}] parameterized by {uij}. For a loop L, the Zi 



flux is given by 4>l = II 



= ±1. 



Defining g = J'/ J, can be diagonalized as 



};s]((A,rV2), (3) 



by finding the complex fermion operators & n £ that are 

linear in Cj's for momentum and band index n = 
1,2,3 (there are six sites per unit cell in the Majorana 
fermion Hamiltonian). This yields the entire spectrum 
£ n feK^y}' s\' ^ ne g roun d states are uniform flux states 
with <fP L — —l for all 12-plaquettes and 4> L = 1 or — 1 for 
all triangular plaquettes, spontaneously breaking time re- 
versal symmetry. A vortex on a plaquette L, defined by 
4>l = —(jP L , costs finite energy for all g. Moreover, the 
uniform flux ground state is degenerate on a torus and 
the topological degeneracy changes across the nA to A 
QPT at <7 C = -\/3CP- However, care is needed for discern- 
ing physical states that satisfy the constraint Eq. (JT|) for 
each configuration of {?% }. 

Topological degeneracy and the projection 
operator- A clue towards an extension of topological 
degeneracy to finite temperature lies in the g-dependent 
effects of the constraint Eq. (111. The constraint de- 
fines the physical states of the free fermion Hamiltonian 
Eq. ([3| and is sensitive to g. Eq. (JT|) can be imple- 
mented using a projection operator P = Yli 1(1 + -D<) 
since D,-P = P [TJ H2] [IB]. Moreover, P commutes with 
the original Hamiltonian 'H[{Uij}]. We can show that 
whether a state survives projection only depends on the 
fermion parity, defined as P/ = Ilijeai'^'z'-imks iCiCj, and 
the parity of the number of vortex excitations |19j on 
triangle plaquettes. 

In the uniform flux sector, all physical states have 
even fermion parity Pf — 1, which is particularly impor- 
tant for determining the topological degeneracy of the 
ground states on a torus. Topological degeneracy comes 
from the identical free fermion spectra, in the thermo- 
dynamic limit, in the four possible topological sectors, 
distinguished by the choice of the Z 2 global flux: 



$ Q EE Yl U \ 



±1, 



(4) 




FIG. 2: A decorated brick wall lattice that is topologically 
equivalent to the star lattice of Ref. pQ, on the surface of a 
torus. Green links denote Uij configurations contributing to a 
global flux threading. The inset defines a links for 'triangles' 
and a' links connecting 'triangles' with a = x,y, z. This la- 
belling of links specifies components of spins interacting across 
the links in the original spin model. 



where a = x,y label two global cycles T x and T y . 
Now ($ x ,<& y ) = (±1,±1) are four distinct states where 
$2, = — 1 and $> y = — 1 indicates tt flux threaded through 
the distinct holes of the torus (see Fig[2]). Normally 
the fermion parity in the unprojected ground state wave 
function in all four topological sectors is even and in- 
dependent of g. Thus, they will survive the projection 
and give rise to four-fold topological degeneracy. This 
is indeed the case on the A side (g > g c ). However, 
the fermion parity in the unprojected ground state wave 
function in the ( — 1,-1) sector on the nA side (<? < g c ) is 
odd and consequently it does not survive projection |20j . 
Thus there is only a three-fold topological ground state 
degeneracy in the nA phase. In summary, for A phases 
with uniform flux, all the physical states consist of an 
even number of fermionic quasiparticle excitations above 
the ground states in all four topological sectors. For nA 
phases with uniform flux, all physical states have an even 
number of fermion excitations in sectors (1,1), (1,-1), 
and (—1,1), but an odd number of fermion excitations 
in the sector (—1,-1). This has consequences not only 
for the topological ground state degeneracy, but also at 
finite temperature, as shown below. 

Global flux expectation value- Motivated by the 
connection between the change in the allowed physical 
spectrum at the topological QPT and the global flux 
states, we consider the 'expectation value' of the global 
flux (<&„) defined as [T2] 



($«(T)) = ~tr $ Q e-«/ T 



(5) 



in a finite size system with N sites. This ties the topo- 
logical degeneracy to the spectrum and offers a natural 
finite T extension of topological degeneracy. If we further 
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restrict ourselves to uniform flux states (which is valid at 
T = and is a good approximation in the vicinity of the 
QCP where the fermion gap vanishes but vortex gap is 
finite), Eq. ^ can be recast as 



E 



* x .*„=±i 



(0) 



where we have defined a sub-partition function for each 
global flux sector (Q x , <& y ). So, for the uniform flux states 



= trl 



(**>*v) 



exp(-W/T). 



Clearly, in the absence of a dependence of the physical 
spectrum on the global flux, all the sub-partition func- 
tions will be identical Z^ 1 ^ 1 ) = = Z^" 1 ) = 
Z^ 1 ' 1 ) and ($ Q ) w iU average out to be identically zero. 
This is the case for the A phase; and the case of toric 
code previously studied [12]. However, for the nA phase 
of the CSL model, the (—1,-1) sector is projected out of 
the ground state Hilbert space and hence Z^~ l ~ x ^{T — 
0) n A = 0. This yields a finite and definite (<fr Q ) in the 
nA phase at T = 0: 



(*x)(T = 0) 



1/3 (nA, 3 < 5c ) 

(A,<?> 5c ) 



(7) 



The significance of Eq.([7]) is that (<5> x ) is the first identi- 
fication of a quantity that can be defined in a thermody- 
namic sense that changes at the topological QPT. Now 
the relation ($> X )(T = 0) can be related to the topological 
degeneracy through 



n DEG = 4-3($ x )(T = 0). 



(8) 



Most importantly, this identification allows one to extend 
the notion of topological degeneracy to finite temperature 
through (<f> x )(T ^ 0) and to investigate the vicinity of the 
topological QPT that is largely unknown in Fig. 1(b). 

Analogous to T c (g) (the solid line of Fig. 1(a)) and the 
cross over line (dashed line in Fig. 1(a)) in the vicinity 
of a conventional QCP, we define and investigate a cross 
over temperature scale T*(g) above which ($> x ) falls off to 
zero for a system with finite size. As shown in Fig. ([3|, a 
crossover temperature scale T*(g) can be defined as the 
temperature scale at which ($> X )(T) falls off exponen- 
tially from its zero temperature value at a given value of 
g < g c . In Fig. pi, we show (<& x ) (T) for g = 1.3 as de- 
fined in Eq.([5]). For all g < g c , (& x ) is nearly a constant 
(1/3) for T < T*(g), but decays exponentially to zero 
at higher temperatures. We have defined T*(g) as the 
point at which Z^ 1 ^^*) / Z {1 ^ {T*) = e" 1 by con- 
vention. Above T* , the distinction between the A and 
the nA phase vanishes. The plot of T*(g) in Fig jl] shows 
that T*(g) is a distinct scale which is non- vanishing at 
the QCP g = g c , unlike the excitation gap which van- 
ishes. The excitation gap is non-vanishing in both phases 
where, in contrast, T*(g) is only non-zero for g < g c 




FIG. 3: Defining T* through the exponential decay of the 
($ X (T)). The plot is for g = 1.3 < g c on a 60 x 40 lattice. 




FIG. 4: The T* (solid line) compared to the single fermion 
excitation gap (dashed line) for a 60x40 lattice in the vicinity 
of g c . Although both Abelian and non-Abelian phases are 
gapped, T* is defined only in the non-Abelian phase. 



which allows the identification of the nA phase at finite 
temperature. 

The crossover scale T* has an intriguing system size 
dependence. For a large system of N sites, 



A(ff) 
IniV ' 



(9) 



where A (g) is the energy gap in the fermion spcctrum(scc 
Figj4j . Hence T* — > as N — > oo, but at a rate slower 
than any other quantity in the system. Therefore, the 
distinction between the A phase and the nA phase, while 
strictly vanishing at finite temperature in the thermo- 
dynamic limit, can be meaningful in some sizable range 
of N. Interestingly, this system size dependence bares 
similarity with the crossover scale for the finite tempera- 
ture topological entanglement entropy of the toric code in 
Ref . [10] . It is also reminiscent of the finite temperature 
behavior of an Ising chain of finite length. An Ising chain 
does not order at any finite temperature in the thermo- 
dynamic limit. However, one can define a finite crossover 
temperature scale in a finite size system of length TV by 
comparing the energy cost 2 J for a domain wall and the 
entropic gain of T log N for N possible choice for the po- 
sition of the domain wall. Nevertheless, the present size 
dependence Eq.pl is rather a consequence of (<& X )(T) be- 
ing bounded from below by tanh Ar (A/2T)/3 (see SOM) 
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FIG. 5: The entanglement entropy change due to a vortex 
pair excitation for 30 x 30 Kagome lattice sites on a torus. 
The change at g c is not sharp for this finite size, though it 
approaches the thermodynamic limit value of In 2 in the nA 
phase. 



and it is unrelated to local fluctuations. Possible con- 
nections underlying these apparent similarities and their 
implications for a "quantum critical region" is an open 
question. 

Towards finite T entanglement entropy While 
the topological entanglement entropy [7J [5] of the ground 
state wave function has been widely adopted as a measure 
of topological order, it is does not distinguish the A phase 
from the n A phase in the present case [17] . The topolog- 
ical entanglement entropy 7 is defined as the universal 
constant term in the entanglement entropy in addition 
to the usual term proportional to the perimeter of the 
boundary: S en t = aL — 7. Further it is known that 7 is 
given by the total quantum dimension of the topological 
phase : 7 = In \/J2 a ^cr Since {d a } = {I, v2, 1} in the 
nA phase while {d a } — {1, 1, 1, 1} in the A phase [T7] . 
the topological QPT does not affect 7 which is equal to 
2 in both phases. However, it is possible a finite tem- 
perature extension of this quantity might offer a possible 
distinction between the two phases, as excitations with 
distinct statistics would contribute. 

An extension of the entanglement entropy to T ^ 
must involve the inclusion of thermal excitations. For 
instance, the extension proposed in Ref.[T^ (an alterna- 
tive definition was proposed in Ref. |15j) retains the basic 
form S cn t = — Tr(— pj± In pX) but uses for the p^, a ther- 
malized reduced density matrix 



E x IT 



(10) 



where I^a^s are energy eigenstates. Clearly Eq.(10l re 



duces to the usual definition at T = when only the 
ground state(s) enter the sum. 

While the zero temperature topological entanglement 
entropies of the A and nA phases are identical, the dif- 
ferent excitations of each phase will generically lead to 
different finite-temperature quantities. As a first pass 



through the problem, we consider the change in the en- 
tanglement entropy in the presence of a pair of vortex 
excitations, one each in the two regions A and B whose 
entanglements are under consideration. Fig. (|5| shows 
the result as a function of g. The additional entropy of 
log 2 of the nA phase reflects the double degeneracy asso- 
ciated with Majorana fermion vortex core states respon- 
sible for the non-Abelian statistics of the vortices. It is 
clear that the characteristics of finite-energy excitations 
are important qualities of the topological phases. 

Closing remarks- We studied the nature of the topo- 
logical QPT between a nA phase and an A phase in 
an exactly solvable model using finite temperature ex- 
tensions of two separate measures of topological order. 
The expectation value of the global flux (<!>)(T) is a fi- 
nite temperature extension of the ground state topolog- 
ical degeneracy which clearly changes at the QPT. We 
found that (4>)(T) retains the T = value for g < g c 
up to a crossover temperature scale which decays loga- 
rithmically with the system size. Whether this type of 
crossover is ubiquitous for topological phases in two spa- 
tial dimension is an open question. As a step towards a 
finite temperature extension of 7, which is independent of 
g at T = 0, we considered the effect of a pair of vortices. 
This indicates the possibility that ~f(T) might distinguish 
the nA phase from the A phase. 
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Note that the original representation in Ref. [T] introduces 

I, 1 gauge fields only on z, z' , and the 'cut' links, which 
corresponds to a particular gauge choice here. 

After gauge fixing, the P expansion of Ref. [16] can be 
factorized: P = [1 + Ui D i\ G 



[1 + -P/EL £ a,v^] G : 



where G = [1+Y\ ^•+E i<J - + ' ■ • ] (note Df = 1). 
The exclusion of the (—1,-1) state from the ground state 
for g < g c is tied to the nA statistics of the vortices. 
A global flux threading $ a is equivalent to the proce- 
dure of (i)creating a vortex pair (ii) transporting one vor- 
tex around the loop r a (iii) annihilating the pair[B]. The 
(—1,-1) sector is equivalent to two vortex loops linked 
to each other, which cannot be undone in the nA phase. 



where A(g) and A(g) + W(g) are respectively the min- 
imum and maximum of the energy spectrume n r[g] that 
depends on g. Since A(g) and W(g) remains finite in the 
thermodynamic limit of N — ► oo Eq. (1 1 4h implies 



lim ($ a 

N-MX) 







(15) 



in the thermodynamic limit. Further using 

[tanh(A/2T)] JV » 1 - 2ATexp(-A/T) at large N 
and low T, we find the system size dependence of the 
crossover temperature: 



A(g) 
IniV ' 



(16) 



THE SUPPORTING ONLINE MATERIAL: THE 
SYSTEM SIZE DEPENDENCE OF T* IN THE 
NON-ABELIAN PHASE. 



Since the vortex gap never closes for all values of g, 
we can ignore the vortex excitation in calculating (<f> x ) 
in the vicinity of the QPT: 



7(1,1) 

J CSL 



r(-l-l) 
J CSL 



07(1,1) , 7 (-l -1) • 

0-6 r or T 



(11) 



J CSL ^ ^CSL 

Since in the nA phase, the fermion occupation number 
parity is odd in the ($ x , Q y ) = (—1,-1) sector and even 
in all the other sectors 



7(1.1) , -l) 

^CSL^^CSL 



= 2exp(-£ G /T) JJ[l±exp(-e n g/T)], 

(12) 



n . k 



where Eq is the ground state energy. With Eq.(12l 



we can recast the global flux thermal expectation value 
($*>CT) as 

n„ctanh(e n ,/2T) 



(*«>(r) 



2 + n„,fctanh(e„ £ /2T)' 



(13) 



This implies 



tanh 



A(g) 
2T 



N 



tanh 



A(g) + W(g) 
2T 



N 



(14) 
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